Ochanomizu University Workshop; June 14, 2001 



The Infrared Boundary of Perturbative QCD 

Victor ELIAS [] 



" Department of Applied Mathematics, The University of Western Ontario, London, 

Ontario N6A 5B7 CANADA 



Abstract 

Evidence is reviewed suggesting that QCD remains a perturbative theory with a (rela- 
tively) small coupling constant down to a distinct infrared boundary on perturbative physics, 
a boundary corresponding to the momentum scale associated with a l3-function pole. 
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1 Supersymmetric Gluodynamics 



Our information about the coupling constants characterising known interactions is 
necessarily perturbative. Many QCD calculations, for example, have been performed 
in the MS scheme in which the (]- function for the couplant x = as{n)/TT is known to 
three subleading orders 

/x^— ^ f5{x) = -Pox' - Pix' - P2x' - Psx' - P,x'... (1) 



Po = 11/4 

= 51/8 - 19n//24 

p2 = 2857/128 - 5033n//1152 + 325n5/3456 

/?3 = 114.23033 - 27.133944^/ + 1.582379nJ + 5.856696 • lO^^nJ (2) 

Thus, we can determine with high accuracy how the QCD couplant evolves with /x 
in perturbatively accessible processes. However, the behaviour of the couplant in the 
infrared region, the region where successive terms in the series (1 ) become comparable 
in magnitude, is not at all clear from the truncated perturbative series. 

There is value in having theories in which the ^-function is exactly known. One 
such theory is N = 1 supersymmetric SU{Nc) Yang-Mills theory, a theory of gluons 
and gluinos often denoted as ^'supersymmetric gluodynamics." The ^-function for 
this theory has been computed to all orders by Novikov, Shifman, Vainshtein, and 
Zakharov (NSVZ) using instanton calculus methods it is also obtainable alge- 
braically by requiring that the supersymmetric anomaly multiplet characterising the 
theory satisfy the Adler-Bardeen theorem ^, 

2 dx _ oNSVZf \ _ _^NcX^ I 

^ dfi^ - ^ ^ ' 4 [l-iV,x/2_ 

X = gV^ir' (3) 

The pole characterising the P- function (3) necessarily implies the existence of an 
infrared boundary for the asymptotically-free phase of the couplant Suppose we 
have some sufficiently small initial value xo [xo < 2/Nc\ occurring at momentum 
scale fiQ. We see that the P-function is negative, implying that x{fj,) grows as fi 
decreases until x = 2/Nc, at which point P{x) becomes singular. Alternatively fi 
can be regarded as a function of x, with initial value /iq occurring at xq and with 
an extremum occurring when x = 2/Nc. This extremum is necessarily a minimum, 
a critical value Hc that represents the lower bound on the domain of x{fi). In short, 
supersymmetric gluodynamics is a field theory characterised by a real perturbative 
interaction coupling constant only when fj, > fi^. Not only is fic the infrared boundary 
of this theory, but this boundary is characterised by a finite value of the couplant: 
x{^c) = 2/N,. 
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The behaviour described above is highly suggestive of the quahtative picture we 
have of the strong interactions, whose perturbative character as a theory of quarks 
and gluons abruptly ends at near-GeV hadronic mass scales. At such mass scales, 
strong interactions are described by some effective-Lagrangian phenomenological the- 
ory distinct from perturbative QCD, a theory in which fundamental fields are no 
longer quarks and gluons, but hadrons. 



2 Non-Supersymmetric Gluodynamics 



For a function underlying an incompletely known perturbative series, such as the QCD 
j3-function (1 ), one of the few available techniques for extracting information about 
the existence of poles or zeros of that function is to examine Fade approximants con- 
structed entirely from the known terms of the series. Such approximants may exhibit 
a positive zero which precedes any positive poles, behaviour indicative of a nonzero 
infrared-stable fixed point. Alternatively, a leading pole, e.g., a positive pole which 
precedes any positive zeros in the approximant, is indicative of a (^-function analogous 
to (2), a function characterised by an infrared boundary to perturbative dynamics in 
its asymptotically free phase. While one particular Fade approximant may exhibit 
a spurious pole (usually denoted as a ''defect pole"), the occurrence of leading poles 
in many different Fade approximants all constructed from the same truncated series 
is strong evidence that such a pole genuinely occurs within the function underlying 
that series /||/. [| 

For (Nc = 3) gluodynamics, this truncated series is given by eqs. (1 ) and (2) with 
Uf = 0. The lowest approximants possessing both zeros and poles that are determined 
by the known and first unknown terms of the series (1 ) are 
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-X 



1 - 2.7996a; - 3.7A75x'^ 
1-5.1178X 

1 - 5.9672X 



1 - 8.2854X+ 11.091x2 



(4) 



(5) 



P 



[3|1] 
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-X 



1 + (2.31818 - 0.0087542/54)x + (8.11648 - 0.020294/34)x2 



+ (41.5383 - 0.071053/54)x^ / [1 - 0.0087542/54x] 



(6) 



[2|2]/N_ 11 2 [1 + (13-4026 - 0.027715/34)x + (-22.9153 + 0.032788/54)x2] 
^ ~ [1 + (11.0844 - 0.027715/34)x + (-56.7275 + 0.097036/34)a;2] ^ ^ 

/Jt^l^l = -^x"^ [1 + (9.56239 - 0.022220/54)3;] / [1 + (7.24421 - 0.022220/34)x 

+ (-24.9099 + 0.005151/54)0;^ + (-42.5901 + 0.060939/54)a;' 



demonstration of how Pade approximant poles for a series that differs only infinitesimally 
from a geometric series are able to reproduce the geometric-series pole is presented in ref. [4] . 
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The Maclaurin expansions of the first two approxiniants reproduce the known terms 
~ of the series (1) when Uf = 0. The Maclaurin expansions of the final three 
approxiniants also reproduce the unknown /?4 term of that series. In (4) and (5), 
the first positive pole [x = 0.195 and x = 0.151, respectively] is seen to precede the 
first positive zero [x = 0.264 and x = 0.168]. Remarkably, this behaviour persists 
in (6), (7) and (8), regardless of the value of the unknown coefficient [3^ in the [3- 
function series (1), as is demonstrated graphically in Figs. 3-5 of ref. In (6), 
for example, a positive pole occurs at x = 114.2//54 provided is positive. This 
pole is always smaller than any positive zeros of the degree-3 numerator polynomial. 
Moreover, when jS^ is negative, eq.(6) exhibits no positive zeros or poles. Eqs. (7) 
and (8), which have higher-degree denominator polynomials, exhibit a positive pole 
for all values of (3^ which is seen to always precede any numerator zeros for that same 
value of P^. Thus, for Fade approximants to the (3- function capable of generating 
both zeros and poles, a zero never precedes a pole for Nc = 3 gluodynamics, as would 
be expected if gluodynamics were to have an infrared stable fixed point. Instead, 
all such approximants point strongly toward the existence of a f3-function pole, as 
is known to characterise supersymmetric gluodynamics in the NSVZ renormalization 
scheme. 

Such behaviour also characterises QCD in the 'tHooft (Nc —>■ oo) limit, a "glu- 
odynamics" for any finite choice of Uf. This is evident from the five approximants 
analogous to (4-8) constructed from the known and leading unknown coefficient in 
the j3-function ^ describing the evolution of the (finite) couplant A = Ncas{fJ')/4:7r 
as Nc ^ oo {as 0): 

/z^^^ = --A^ - -A^ - ^A^ - 315.49A^ - /SA'- (9) 

These five approximants are presented in the Appendix to ref. /|^. The pole/zero 
structure of these approximants is the same as that of Nc = 3 gluodynamics - a 
positive pole always precedes any positive zeros occurring within leading approximant 
versions of the f3 -function (9). These results strongly point to the existence of an 
infrared boundary to gluodynamics as a perturbative theory, a boundary which occurs 
at the value of /i associated with such poles. 

3 QCD 

In ref. the dynamics described above for gluodynamics are shown to apply to 
all approximant versions of the QCD (3-function, even if n j is as large as 5. When 
nf >Q, /^[^l^l no longer exhibits a pole, and poles are seen to precede zeros in 
/S^^'^l and /^[^'^l only if (3/^ is larger than approximant-specific lower bounds. Such 
results corroborate a lattice study indicative of a similar Uf threshold {uf = 7) for 
infrared-stable fixed points to occur within QCD 

Fade approximant methods have been used successfully to predict unknown QCD 
(3-function terms. For example, such methods, accompanied by explicit knowledge of 
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the 0{n^j) term in led Ellis, Karliner and Samuel ^ to predict in (1): 

, _ 23,600- 6400n^ + 350n^ + 1.499n3 

~ 256 • ^^"^ 

Comparison of this result with the exact result (2), as calculated explicitly in 
demonstrates the power of Fade approximant methods. Similar methods have been 
used to obtain a corresponding prediction for the presently-unknown coefficient 13^ 

p^grf _ 759, 000 - 219, OOOn/ + 20, 500n^ - 49.8^1^ - 1.84n) 
= 1024 • 

In attempting to extract information about the infrared boundary of perturbative 
strong interaction physics, our interest is necessarily directed toward the Uf = 3 case 
of (11), for which (3^^'^ = 278. The Fade approximant versions of the nj = 3 QCD 
MS (3-function whose Maclaurin expansions reproduce this estimate of [3^ as well as 
j3i — Ps, as given in (2), are 

^r.ni, , 9 2[l-4.116x- 6.006a;2- 5.359a;31 

= — x'^- (12) 

^ ^ ^ 4 1 - 5.893X ^ ^ 

p,,] _ 9 [l- 5.498x- 1.972x^] 
^ ^ ' ~ 4 [1 - 7.276a; + 6.492x2] ^ 

^ ^ ^ 4 [1 - 7.517X + 8.893x2 -3.190x3] v ; 

All three approximants exhibit positive poles that precede any positive zeros. These 
poles Xc respective occur at 



X 



f'^^ = 0.170, xt^l^] = 0.160, x'^l^l = 0.162 (15) 



and correspond to values of between 0.50 and 0.53. 

Both the agreement exhibited by these distinct approximants on the magnitude 
of their poles, as well as the perturbatively small value of the a.^ corresponding to 
these poles, are quite striking. The mass scale fic associated with these approximant 
poles is easily obtained from integrating the differential equation 

/^^^x = /3[^l*^](x) (16) 
from a physical initial value Xo(/io) to the pole value xj.^'*^': 

(17) 



1 dx 



2 7.0 /^i^i^'^ii 



X 



If we identify as{m^) = 0.314 ± 0.010 jl^/ and if we require the couplant to devolve 
with four active flavours from /i = m,- to a four- flavour threshold at 1.2 GeV (below 
which only three active flavours contribute), we find that 

^[3|i] = 950 ± 50 MeV, /i^^l^] = 990 ± 50 MeV, /i^^] = 980 ± 50 MeV. (18) 
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Alternatively, we can utilise the (three-active flavour) estimate asimr) = 0.33 ± 0.02 
to Gnd entirely via (12-14) and (17) that g: 



/ii^l^l = 1.09 ± 0.11 Gel/, /i^l^l = 1.14 ± 0.11 GeV, /i^^^ = 1.13 ± 0.11 GeV. (19) 

These estimates are all suggestive of the fic — 47r/^ (=1.17 GeV) Georgi-Manohar 
boundary on effective-Lagrangian strong interaction physics These results also 
show remarkable consistency across the approximants (12-14) in the coordinates of 
the infrared terminus (fj^c^c = c^sifJ'c)/'^) of QGD as a perturbative theory of quarks 
and gluons. 

Of course, these estimates rely on the estimate (11) for (3^. Because of the alter- 
nation of sign in this expression, the estimate itself is considerably less likely to be 
accurate than the individual estimates for polynomial coefficients of powers ofuf. An 
alternative approach toward acquiring insight into the infrared boundary of QGD is 
to assume that as{fic) = 7r/4, the threshold value associated with dynamical chiral 
symmetry breaking Individual [3|1], [2|2], and [1|3] approximants to the Uf = 3 
(3-function when (3^ is taken to be arbitrary are given by 
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= -^x^ [i + (1.7778 - 0.021 174/54)a; + (4.4471 - 0.037642/34)a; 

+ (20.990 - 0.094670/?4)a;^] / [1 - 0.021174/54x] (20) 

[2|2]/x_ 9 + (7.19456 -0.045604/34)x + (-11.3292 + 0.033619/54)x2] 
P [x) - [1 + (5.41678 - 0.045604/34)a; + (-25.4301 + 0.11469/34)a;2] ^ ^ 

/3[i|3l(a;) = [1 + (5.80845 - 0.041491/54)x] / [1 + (4.03067 - 0.041491/34)a; 

+ (-11.6367 + 0.073762/34)x2 + (-18.3242 + 0.054377/54)x^] (22) 

These approximants are respectively seen to develop a leading positive pole atx = 1/4 
(i.e. as = n/4:) for the following values of (3^: 

/3f '^1 = 189, /3f = 181, /3f '^^ = 202, (23) 

values surprisingly uniform across the three distinct approximants. Upon substituting 
such values for (3^ into their corresponding approximants (20-22), we can estimate via 
(17) concomitant values for the infrared-boundary momentum scale at the x = 1/4 
pole. We set x\!^^^'^^ = 1/Ain (17) and utilize the extreme values of the 3-active-fiavour 
range 0.33 ± 0.02 already quoted for asim^-) in order to find that [11] 



a,(m,) = 0.31 : = 775 MeV, /if = 785 MeV, /i^^' = 778 MeV 
as{mr) = 0.33 : fif^^^ = 863 MeV, /if = 874 MeV, jj.™ = 867 MeV 
as{mr) = 0.35 : /tf = 948 MeV, /tf ^1 = 960 MeV, /tf^l = 952 MeV 
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There is remarkable uniformity in these estimates obtained from distinct approxi- 
mants, suggesting that the infrared boundary is genuine, rather than a single- approximant 
artefact. A more precise estimate of the boundary, however, requires more precise 
knowledge of ag at fi = m^. We can conclude, however, that there is clear evidence 
from Pade approximant methods for QCD to be a computationally-perturbative gauge 
theory of quarks and gluons right up to an 0{1 GeV) boundary mass scale, an ex- 
plicit lower bound on the domain of as{iJ>) below which the description of the strong 
interactions must necessarily be quite different. 
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